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The current status and prospects of research on ferromagnetism in semiconductors are reviewed. 
The question of the origin of ferromagnetism in europium chalcogenides, chromium spinels and, 
particularly, in diluted magnetic semiconductors is addressed. The nature of electronic states derived 
from 3d of magnetic impurities is discussed in some details. Results of a quantitative comparison 
between experimental and theoretical results, notably for Mn-based III-V and II- VI compounds, 
are presented. This comparison demonstrates that the current theory of the exchange interactions 
mediated by holes in the valence band describes correctly the values of Curie temperatures Tq, 
magnetic anisotropy, domain structure, and magnetic circular dichroism. On this basis, chemical 
trends are examined and show to lead to the prediction of semiconductor systems with Tq that may 
exceed room temperature, an expectation that are being confirmed by recent findings. Results for 
materials containing magnetic ions other than Mn are also presented emphasizing that the double 
exchange involving hoping through d states may operate in those systems. 



INTRODUCTION 

Because of complementary properties of semiconduc- 
tor and ferromagnetic material systems, a growing ef- 
fort is directed toward studies of semiconductor-magnetic 
nanostructures. Applications in in optical modulators 
and insulators, in sensors and memories §, as well as 
for computing using electron spins can be envisaged. 
The hybrid nanostructures, in which both electric and 
magnetic field are spatially modulated, are usually fab- 
ricated by patterning of a ferromagnetic metal on the 
top of a semiconductor || or by inserting ferromagnetic 
dots or layers into a semiconductor matrix In such 
devices, the stray fields can control charge and spin dy- 
namics in the semiconductor. At the same time, spin- 
polarized electrons in the metal can be injected into or 
across the semiconductor Furthermore, the ferro- 
magnetic neighbors may affect semiconductor electronic 
states by the ferromagnetic proximity effect even under 
thermal equilibrium conditions. Particularly perspective 
materials in the context of hybrid structures appear to 
be those elemental or compound ferromagnets which can 
be grown in the same reactor as the semiconductor coun- 
terpart. 

However, already the early studies of Cr spinels ||, 
rock-salt Eu- @, f| |, || and Mn-based jl| chalco- 
genides led to the observation of a number of outstand- 
ing phenomena associated with the interplay between 
ferromagnetic cooperative phenomena and semiconduct- 
ing properties. The discovery of ferromagnetism in zinc- 
blende III-V U|, g and II- VI H, |l§ Mn-based com- 
pounds allows one to explore physics of previously not 
available combinations of quantum structures and mag- 
netism in semiconductors. For instance, a possibility 
of changing the magnetic phase isothcrmally, by light 
(l|, [l?], Q or by the electric field (If, |(1 was put 
into the evidence in (In,Mn)As/(Al,Ga)Sb |L7|, and 
(Cd,Mn)Te/(Cd,Zn,Mg)Te ||, |l§ |o) heterostructures. 



The injection of spin-polarized carriers from (Ga,Mn)As 
to a (In,Ga)As quantum well in the absence of an external 
magnetic field was demonstrated, too At the same 
time, outstanding phenomena, known from the earlier 
studies of metallic multilayer structures, have also been 
observed in ferromagnetic semiconductors, examples be- 
ing interlayer coupling [^2[ [||, exchange bias [2J|, giant 
|p2| and tunneling magnetoresistance [^5| . It is then the 
important challenge of materials science to understand 
the ferromagnetism in these compounds and to develop 
functional semiconductor systems with the Curie temper- 
atures Tq exceeding comfortably the room temperature. 

We begin this review by describing briefly, in Sec. 2, 
various families of ferromagnetic semiconductors and the- 
oretical models proposed to explain the nature of rele- 
vant spin-spin exchange interactions. We then discuss, 
in Sec. 3, energies and character of electronic states de- 
rived from 3d shells of magnetic impurities in II- VI and 
III-V semiconductors, which provide information on the 
charge and spin states as well as on the effect of the mag- 
netic constituent on the carrier concentration, also in the 
presence of co-doping by shallow donors or acceptors. In 
Sec. 4, we outline the main ingredients and limitations of 
the mean-field Zener model put recently forward to de- 
scribe quantitatively the hole-mediated ferromagnetism 
in tetrahedrally coordinated semiconductors ^(| ^8) . 
Results of a quantitative comparison between experimen- 
tal and theoretical results for Mn-based III-V and II- VI 
compounds are shown in Sec. 5. This comparison demon- 
strates that the current theory of the exchange interac- 
tions mediated by holes in the valence band describes cor- 
rectly the values of Tq, magnetic anisotropy, and domain 
structure. Finally, in Sec. 6, we present the theoretically 
predicted chemical trends and discuss various suggestions 
concerning the design of high-temperature ferromagnetic 
semiconductors . 

The novel physics of magnetic heterostructures consti- 
tutes the topic of the next paper |29|| in this compendium 
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and, therefore, is not discussed here. There are, of course, 
a number of other recent review articles describing many 
aspects of ferromagnetic III-V @ |l[ ||, II- VI || ||, 
and IV- VI materials [[l2| , which are not touched upon in 
this paper. 



FAMILIES OF FERROMAGNETIC 
SEMICONDUCTORS AND RELEVANT 
SPIN-SPIN EXCHANGE INTERACTIONS 



Manganites (perovskite: (La,Sr)MnC>3 and related ma- 
terials), which show colossal magnetoresistance, are mag- 
netic semiconductors, whose studies have been particu- 
larly active over the recent years |35[ p6fl . Their ferro- 
magnetic order, beginning at 350 K, is brought by the 
double-exchange interaction involving on-site Hund's fer- 
romagnetic spin coupling and hopping of d electrons be- 
tween Mn 3+ and Mn 4+ ions. 

The family of magnetic semiconductors encompasses 
also europium and chromium chalcogenides (rock-salt 
type: EuS, EuO and spinels: CdCr 2 S 4 , CdCr 2 Se 4 ), for 
which the Curie temperature Tq does not exceed 100 K. 
In the case of rock-salt Eu compounds, there appears 
to be a competition between antiferromagnetic cation- 
anion-cation and ferromagnetic cation-cation superex- 
change 0|. The latter can be traced back to the ferro- 
magnetic s-f coupling, and the presence of s-f hybridiza- 
tion, which is actually stronger than the p-f hybridization 
due to symmetry reasons j9| |37j. In such a situation, 
the lowering of the conduction band associated with the 
ferromagnetic order enhances the energy gain due to hy- 
bridization. The Cr-spinels represents the case, in which 
the d orbitals of the two cations are not coupled to the 
same p orbital, which results in a ferromagnetic superex- 
change. 

In europium and chromium chalcogenides discussed 
above, the presence of carriers can affect Tq but is not 
necessary for the appearance of the ferromagnetic order. 
In contrast, in the case of Mn-based IV- VI Jl|, III-V [f|§, 
and II- VI ||{| diluted magnetic semiconductors (DMS) 
the ferromagnetism can be observed provided that the 
hole concentration is sufficiently high. According to ex- 
perimental and theoretical results, which are summarized 
in next sections, the ferromagnetic order in Mn-based 
DMS is mediated by carriers residing in relatively wide 
valence bands. In the case of III-V and II- VI DMS, the 
holes are coupled to the localised spins via a strong, sym- 
metry allowed, antiferromagnetic p-d interaction. This is 
in contrast to manganites, where only d electrons in nar- 
row bands appear to be involved. However, the d levels of 
transition metals other than Mn reside in the band gap of 
III-V and II- VI compounds. In such a situation the dou- 
ble exchange may constitute the dominant mechanism of 
spin-spin interactions. 
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FIG. 1: Approximate positions of transition metals levels rel- 
ative to the conduction and valence band edges of II- VI (left 
panel) and III-V (right panel) compounds. By triangles the 
d^/d^ -1 donor and by squares the d N /d N+1 acceptor states 
are denoted (adapted from Ref. [|43[). 



MAGNETIC IMPURITIES IN II- VI AND III-V 
COMPOUNDS 

Energy levels 

We consider tetrahedrally coordinated semiconductors, 
in which the magnetic ion occupies the cation sublat- 
tice, as found by extended x-ray absorption fine structure 
(EXAFS) studies in the case of Cdi_ x Mn x Te go| and 
Gai-^Mn^As Obviously, magnetic properties of a 

semiconductor containing magnetic ion will depend on 
energetic positions of states derived from magnetic shells 
as well as on their interactions with the host bands. Fur- 
thermore, the energy of the magnetic levels in respect to 
host bands together with the on-site correlation energy U 
determine whether the magnetic ion act as a dopant and 
how its charge and spin state depends on the presence of 
other impurities. 

According to the internal reference rule p^ |, the posi- 
tions of states derived from magnetic shells do not vary 
across the entire family of the II- VI or III-V compounds if 
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the valence band offsets between different compounds are 
taken into account. In Fig. 1, adapted from Ref. jl3|], the 
data for II- VI and III-V DMS containing various tran- 
sition metals are collected. The symbols D(0/+) and 
A(0/— ) denote the donor and acceptor levels derived 
from 3d shells of magnetic ions. In the case of II- VI DMS 
(left panel), these states correspond to the transforma- 
tion of the doubly ionised magnetic ions M 2+ into M 3+ 
and into M 1+ ions, in their ground states, respectively, 
that is to the lower and upper Hubbard bands. Similarly, 
in the case of III-V DMS (right panel) D(0/+) and A(0/- 
) denote the donor and acceptor states which, however, 
in contrast to the situation in II- VI DMS, correspond 
to the transformation of the triply ionised magnetic ions 
M 3+ into M 4+ and into M 2+ ions, respectively. A charac- 
teristic evolution of the level energies with the number of 
the d electrons is seen in both II- VI and III-V DMS, the 
pattern known from atomic spectra but significantly flat- 
tened out in solids by screening and hybridization effects 
S 

It should be noted at this point that the internal refer- 
ence rule may serve only for the illustration of chemical 
trends and not for extracting the precise values of the ion- 
isation energies. Moreover, the interaction between the 
impurity and host states can lead to the appearance of 
additional band-gap levels derived from the semiconduc- 
tor bands. This may cause some ambiguity concerning 
the nature of localised states observed experimentally in 
these systems. In particular, a strong p-d hybridization 
can lead to a binding of a hole in a Zhang-Rice polaron 
(charge transfer) state, which then gives rise to an ad- 
ditional level in the band gap ||], |(| |47[ ||, ||| . Fur- 
thermore, if the d^/d^ -1 donor state resides above the 
bottom of the conduction band, the ground state corre- 
sponds to a hydrogenic-like level d 7V_1 +e located below 
the band edge, as observed in CdSe:Sc Q. Similarly, if 
the acceptor state lies under the top of the valence band, 
the ground state corresponds to a hydrogenic-like accep- 
tor d w+1 +h, not to the d N state. Importantly, band 
carriers introduce by such magnetic ions can mediate ex- 
change interactions between the parent spins. Obviously, 
energies of hydrogenic-like states follow the band edges, 
and by no means are described by the internal reference 
rule. This appears to be the situation of the Mn related 
levels in the gap of III-V compounds (IsA , the case dis- 
cussed in detail below. 



Mn in II- VI compounds 

It is well established that Mn is divalent in II- VI com- 
pounds, and assumes the high spin d 5 configuration char- 
acterized by S = 5/2 and g = 2.0 @ §| ||. Indeed, ac- 
cording to Fig. 1, the Mn ions neither introduce nor bind 
carriers, but gives rise to the presence of the localised 
spin in II- VI DMS. The spin dependent hybridization 



between anion p and Mn d states leads to the superex- 
change, a short-range antiferromagnetic coupling among 
the Mn moments. In order to take the influence of this 
interaction into account, it is convenient to parameterize 
the dependence of magnetization on the magnetic field 
in the absence of the carriers, M a (H), by the Brillouin 
function, in which two empirical parameters, the effec- 
tive spin concentration x e qNo < xNq and temperature 
T e jj > T, take the presence of the superexchange interac- 
tions into account |52| [53| . The dependencies x e g(x) and 
T^p(x) = T e ff(x) — T have been determined for a number 
of Mn-based II- VI DMS. Importantly, the antiferromag- 
netic superexchange can be overcompensated by ferro- 
magnetic interactions mediated by band holes (2(|, the 
theoretical prediction confirmed subsequently by the ob- 
servation of ferromagnetic ordering in p-type II- VI DMS 

MM. 



Mn in III-V compounds 

Figure 2, taken from Ref. (49), shows the energetic po- 
sition of the Mn impurity level in III-V compounds, as 
evaluated by various authors from measurements of opti- 
cal spectra and activation energy of conductivity. A pri- 
ori, the Mn atom, when substituting a trivalent metal, 
may assume either of two configurations: (i) d 4 or (ii) 
d 5 plus a weakly bound hole, d 5 +h. Accordingly, the 
experimentally determined energies correspond to either 
d 4 /d 5 or d 5 +h/d 5 levels. 

It appears to be a general consensus that the Mn 
acts as an effective mass acceptor (d 5 +h)in the case of 
antimonides and arsenides. Such a view is supported 
by the relatively small Mn concentrations leading to 
the insulator-to-metal transition, which according to the 
Mott criterion n 1 / 3 as = 0.26, points to a relatively large 
extension of the effective Bohr radius as- Moreover, 
the ESR studies of GaAs:Mn reveal, in addition to the 
well known spectrum of Mn d 5 with the Lande factor 
9Mn = 2.0, two additional lines corresponding to g\ = 2.8 
and (72 ~ 6 |54|, |55|, [S6| , which can be described quantita- 
tively within the k ■ p scheme for the occupied acceptor 
f34] , H). Here, the presence of a strong antiferromag- 
netic p-d exchange interaction between the bound hole 
and the Mn d-electrons has to be assumed, so that the 
total momentum of the complex is J = 1. In agreement 
with the model, the additional ESR lines, in contrast 
to the g — 2.0 resonance, are visible only in a narrow 
range of the Mn concentration |36| . This range should be 
greater than the concentration of compensating donors, 
and smaller than that at which acceptor wave functions 
start to overlap and merge with the valence band. The 
antiferromagnetic coupling within the d 5 +h complex is 
seen in a number of experiments, and has been employed 
to evaluate the p-d exchange integral f3N ~ -1 eV Q 
in GaAs:Mn, the value in agreement with that deter- 
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FIG. 2: Experimental energies of Mn levels in the gap of 
III-V compounds according to in respect to valence-band 



edges,whose relative positions are taken from Ref. [151 1 (after 
Ref. @. 



mined from interband magnetoabsorption in (Ga,Mn)As 

Importantly, the above scenario is corroborated by re- 
sults of photoemission [ |59| , |60| and x-ray magnetic cir- 
cular dichroism (XMCD) studies |il], in metallic or 
nearly metallic (Ga,Mn)As. The latter point to the d 5 
Mn configuration. The former are not only consistent 
with such a configuration but also lead to the value of 
(3Nq similar to that quoted above, (3Nq w —1.2 eV [ [d9| . 
Furthermore, the photoemission reveals the presence of 
two features in the density of states brought about by 
the Mn constituent: the original Mn 3d 5 states located 
around 4.5 eV below the Fermi energy £/p, and new states 
merging with the valence band in the vicinity of Ep |io) . 
These new states correspond to acceptors, as discussed 
above. They are derived from the valence band by the 
Coulomb field as well as by a local Mn potential that 
leads to a chemical shift in the standard impurity lan- 
guage, or to a valence band offset in the alloy nomencla- 
ture. 

In contrast to antimonides and arsenides, the situation 
is much more intricate in the case of phosphides and ni- 
trides. Here, ESR measurements reveal the presence of 
a line with g = 2.0 only |63, 64, 
assig ncd to d 5 centers @ |6j, |66| 



5| £6|, which is thus 
Moreover, according 



to a detailed study carried out for a compensated n-type 
GaP:Mn j6f|, the ESR amplitude diminishes under illu- 
mination and, simultaneously, new lines appear, a part of 
which exhibit anisotropy consistent with the d 4 configu- 
ration. This, together with the apparent lack of evidence 
for d 5 +h states, even in p-type materials, seems to im- 
ply that Mn in the ground state possesses the d 4 , not 
d 5 +h, electron configuration |6f|. This would mean that 
the Mn energy in Fig. 2 for GaP fiBl and, therefore, for 



GaN |67], ^8) (where the valence band is lower than in 
GaP) corresponds to the d 4 /d 5 , not d 5 +h/d 5 level. Such 
a view appears to be supported by the ab initio computa- 
tion within the local spin density approximation (LSDA) , 
which points to the presence of the d-states in the gap 
of (Ga,Mn)N |39). In this situation, the spin-spin inter- 
action would be driven by a double exchange mechanism 
involving hopping of d-electrons |6^, [7fJ , as in the case of 
colossal magnetoresistance manganites, not by the holes 
in the valence band. 

However, the above interpretation has recently been 
call into question In particular guided by photoe- 
mission results for II- VI compounds |4q| one expects that 
the energy of the d 4 /d 5 level will vary little between ar- 
senides and nitrides. This implies that this level should 
reside in the valence band of GaN despite the 1.8 eV va- 
lence band offset between GaN and GaAs, as shown in 
Fig. 1. The resulting contradiction with the LSDA find- 
ings can be removed by noting that in the case of strongly 
correlated 3d electrons, a semi-empirical LSDA+U ap- 
proach is necessary to reconcile the computed and pho- 
toemission positions of states derived from the Mn 3d 
shell in (Ga,Mn)As @ 0. 

Another important aspect of magnetic acceptors is 
that the p-d hybridization, in addition to producing the 
exchange interaction, can contribute to the hole binding 
energy Eb- By taking the hole wave function as a coher- 
ent superposition of p-states of anions adjacent to Mn |4[| 
and assuming the p-orbitals to be directed towards the 
Mn ion, the state has been found to have 30% lower 
energy than that corresponding to the mutually paral- 
lel p-orbitals E{|. This shows that the Kohn-Luttinger 
amplitudes away from the T point of the Brillouin zone 
are also involved. In order to evaluated Eb, a square- well 
spherical potential V(r) = V Q(b — r) is assumed J17), 
whose depth V Q is determined by the p-d hybridization 
taking into account the above mentioned arrangement of 
the p-orbitals, 



V a = 



5 PN 
8tt 1.04 



1 - 
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Here, the values of (3 No are taken from Ref. p8| ; A e g- 
is the distance of d 4 /d 5 level to the top of the valence 
band, which is evaluated to be 2.7 eV in (Ga,Mn)As fl59| , 
and is assumed to be reduced in other compounds by the 
corresponding valence band offsets, and = 7 eV is 
the correlation energy of the 3d electrons |59|, the 
energies visible in Figs. 1 and 2. Finally, b/a Q is the well 
radius in the units of the lattice constant, and should lie 
in-between the cation-anion and cation-cation distance, 
Voa /4 < b < a /v2. It turns out that in the case of 
GaN:Mn the hole is bound by Mn, even in the absence 
of the Coulomb potential, Eb = 1.0 eV for b — 0.46a o . 
This demonstrates rather convincingly that a large part 
of Eb originates indeed from the p-d interaction, indicat- 
ing that the Zhang-Rice (ZR) limit Esl [47| is reached in 
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these systems. One of the important consequences of the 
the ZR polaron formation is the shift of the Mott critical 
concentration towards rather high values. According to 
the known relation between and as |72| , the critical 
hole concentration is p c = 4 x 10 19 cm^Mn (Ga,Mn)As 
and at least an order of magnitude greater in (Ga,Mn)N, 
if no shallower acceptors are present. 



THEORETICAL DESCRIPTION OF MN-BASED 
FERROMAGNETIC SEMICONDUCTORS 

Electronic states and models of carrier-mediated 
ferromagnetism 

Despite a considerable effort aiming at elucidating 
the nature of ferromagnetism in Mn-based II- VI and 
1II-V DMS, the form of the relevant minimal Hamil- 
tonian and its universality for all compounds are still 
under an active debate. Such a situation reflects the 
multifaceted environment, in which the ferromagnetism 
appears. Indeed, conceptual and technical difficulties 
inherent to theory of strongly correlated and disor- 
dered transition-metal compounds are combined — in fer- 
romagnetic semiconductors — with the intricate physics of 
Anderson-Mott localisation and disordered Stoner mag- 
netism that are specific to heavily doped semiconductors. 
Moreover, low-temperature epitaxy, by which III-V ma- 
terials are obtained, results in a large concentration of 
native defects such as antisites, which act as compensat- 
ing donors. Another possible source of local bonds recon- 
structions is the mechanism of self-compensation, occur- 
ring in heavily doped semiconductors once the Fermi level 
reaches the energy triggering defect reactions. Structural 
faults may form with neighbor transition metal impu- 
rities defect complexes exhibiting hitherto non-explored 
magnetic characteristics. At the same time, strong com- 
pensation by donor-like defects enhances the electrostatic 
disorder substantially, leading to deep and long-range po- 
tential fluctuations that result in significant band tailing. 

There are two basic approaches to theoretical mod- 
eling of the materials in question: (i) ab initio or first 
principles studies and (ii) theories starting from effec- 
tive Hamiltonians containing experimentally determined 
parameters. Since ab initio works are reviewed else- 
where j73|, they will not be discussed here except for 
noting that they raise interesting questions to what ex- 
tend the local spin density approximation (LSDA) can 
handle strong correlation inherent to charge-transfer and 
Hubbard-Mott insulators as well as whether the coher- 
ent potential approximation (CPA) can capture the key 
aspects of the Anderson-Mott localisation. Our main fo- 
cus will be on presenting models employing parametrized 
Hamiltonians, and on showing their capabilities and lim- 
itations vis a vis experimental results. 

As discussed above, various experiments imply consis- 



tently that the Mn ions are in the high spin 2+ charge 
state in tetrahedrally coordinated magnetic semiconduc- 
tors. Accordingly, the Mn ions are electrically neutral in 

II- VI compounds but act as effective mass acceptors in 

III- V semiconductors. It is now well established that in 
the absence of free carriers the dominant exchange mech- 
anism is the short-range superexchange in zinc-blende 
magnetic semiconductors. This mechanism leads to anti- 
ferromagnetic interactions, except perhaps for some Cr- 
and V-based compounds, for which the presence of a fer- 
romagnetic coupling is theoretically predicted |74, [75] . 
Remarkably, owing to the large exchange energy (3Nq\ 
and the high density of states, the hole-mediated long- 
range ferromagnetic exchange interaction can overcome 
antiferromagnetic superexchange Indeed, as al- 
ready emphasized, the presence of holes is essential for 
the existence of the ferromagnetic order in Mn-based 
semiconductors. Furthermore, the relevant spin-spin in- 
teraction is long range according to neutron studies |76) . 

It should be recalled at this point that electronic states 
in doped semiconductors undergo dramatic changes as a 
function of the impurity concentration ]77], f78|. Hence, 
the hole states, and possibly hole-mediated exchange 
mechanisms, may a priori undergo dramatic changes as 
function of the Mn content x and the concentrations of 
acceptors Na and compensating donors Nrj. The evolu- 
tion of electronic states in doped semiconductors is gov- 
erned by the ratio of the average distance between the 
carriers r c to the effective impurity Bohr radius as, de- 
termined by both Coulomb field and short-range poten- 
tial of Eq. 1. In the case of the holes in (Ga,Mn)As, 
r c = (3/471-p) 1 / 3 , p = xN - N D , and a B « 0.78 nm 
| f72[ . A similar value is expected for (Zn,Mn)Te contain- 
ing nitrogen acceptors p(jj |. However, as already men- 
tioned, the p-d interaction may significantly contribute 
to the impurity binding energy and diminish the effec- 
tive Bohr radius. In the range of small impurity con- 
centrations, r c 3> as, the holes are tightly bound to ac- 
ceptors. Hence, the conductivity vanishes in the limit of 
zero temperature. At non-zero temperatures, the charge 
transport proceeds either via phonon-assisted hopping 
between occupied and empty acceptors or by means of 
thermal activation from the acceptor levels to the valence 
band. In a pioneering work Pashitskii and Ryabchenko 
|79| evaluated the strength of exchange interactions be- 
tween localised spins mediated by band carriers thermally 
activated from impurity levels. More recently, Wolff et 
al. [|80t ] considered carriers localised on impurities and 
forming bound magnetic polarons (BMP). It was found 
that there exists a range of parameters, in which the cou- 
pling between the BMP is ferromagnetic. This idea was 
further explored by Bhatt and Wan |8lJ], who examined 
by Monte Carlo simulations properties of a ferromagnetic 
phase transition driven by the interactions between the 
BMP. In a more recent work Berciu and Bhatt |7^] dis- 
cuss, within the MFA, ferromagnetic interactions medi- 
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ated by quantum hopping of holes within the Mn accep- 
tor impurity band in III-V materials. The compensation 
is taken into account by assuming that the number of 
the holes is smaller than that of the Mn sites, but no 
effects of ionised donors on the site energies is taken into 
account. An important result is that the disorder in im- 
purity positions tends to enhance the magnitude of Tq. 

Two other groups noted that a long-range exchange 
interaction between Mn spins can be mediated by holes 
undergoing quantum hoping from the Mn-derived im- 
purity states to the extended valence band states. In- 
oue et al. Q adopted the Slater-Koster approach, well 
known in the physics of resonant states, for the case of 
two magnetic impurities. It has been found, by a model 
calculation, that the pairs of Mn spins coupled to the 
valence band states have a lower energy in the ferromag- 
netic than in the antiferromagnetic configuration. Litvi- 
nov and Dugayev ||| suggested than the ferromagnetic 
spin-spin interaction can originate from virtual excita- 
tions between the acceptor-like impurity level and the 
valence band, a variant of the Bloembergen- Rowland in- 
direct exchange mechanism. They evaluated Curie tem- 
peratures by using a formula, derived originally for exci- 
tations between valence and conduction bands, without 
proving its correctness for the case in question. 

With the increase of the net acceptor concentration, 
the impurity band merges with the valence band. For 
r c <C as, the holes reside in the band, and their quasi- free 
propagation is only occasionally perturbed by scatter- 
ing of Mn and other defect potentials, whose long-range 
Coulomb part is screened by the carrier liquid. Here, 
the celebrated Ruderman-Kittel-Kasuya-Yosida (RKKY) 
mechanism, driven by intraband virtual excitations, is 
expected to dominate. In the context of III-V magnetic 
semiconductors, this mechanism was discussed by Gum- 
mich and da Cunha Lima [Q and Matsukura et al. [§51. 
At the same time, the present author and co-workers pq 
demonstrated the equivalence of the RKKY and Zener 
p6| , |87| models, at least on the level of the mean-field 
and continuous medium approximations. However, with 
no doubts, beyond those approximations such equiva- 
lence can be questioned |Q . Within the Zener approach 
, [S7|] , and its nuclear spin variant J89|, the degree of 
spin ordering, M q , at given temperature T is sought by 
minimizing the total free energy of the spin and carrier 
subsystems, F[M q }. Here, M q denotes the Fourier com- 
ponents of localised spin magnetization M(r), so that 
the minimum of F[M q ] for M q= o 7^ implies the fer- 
romagnetic order. Because of its relevance, the Zener 
model will be discussed in some details in a subsequent 
subsection. 

In view of the above discussion a question arises 
whether the hole-mediated ferromagnetism appear in the 
insulator or in the metallic phase. It is well estab- 
lished that the metal-insulator transition (MIT) occurs 
at r c 2Aas in doped non-magnetic semiconductors 



. According to this criterion one gets the critical hole 
concentration p c = 4 x 10 19 cm -3 for as — 0.78 nm. 
Experimentally, the MIT occurs at about 3.5% of Mn in 
(Ga,Mn)As, i.e., for N x = 7 x 10 20 cm" 3 @ ||. A 
large difference between these two values is presumably 
caused by the compensation as well as by the enhance- 
ment of localisation by the sp-d exchange scattering [^3) , 
an effect observed also in p-(Zn,Mn)Te llfj]. This is docu- 
mented in both (Ga,Mn)As @ and p-(Zn,Mn)Te Jl§ by 
the presence of negative magnetoresistance and the asso- 
ciated insulator-to-metal transition driven by the mag- 
netic field In addition to the MIT ati« 0.035 in 
(Ga,Mn)As, a reentrant insulator phase is observed for 
x > 0.06 65). Presumably, a self-compensation mech- 
anism is involved, such as the appearance of interstitial 
Mn donors, as suggested by first principles studies Q. 
Further efforts are necessary to test this hypothesis and, 
more importantly, to push the Mn solubility limit to- 
wards higher x values. 

Perhaps, the most intriguing property of the materials 
in question is that the ferromagnetism is observed on the 
both sides of MIT |9l{| . It is, therefore, interesting 

to contemplate the nature of electronic states in the vicin- 
ity of the MIT in doped semiconductors. Obviously, the 
random spatial distribution of acceptor and donor cen- 
ters gives rise to strong spatial fluctuations in the carrier 
density and states characteristics. According to the phe- 
nomenological two-fluid model there exist two kinds of 
relevant states |)4|. The first are strongly localised and 
thus singly occupied states associated with the attractive 
potential of a single majority impurity. The strongly lo- 
calised carriers barely contribute to the conduction pro- 
cess. However, they produce a Curie-like component in 
the magnetic susceptibility and give rise to the pres- 
ence of BMP in magnetic semiconductors. Obviously, 
the impurity-like states dominate deeply in the insulating 
phase but their presence is noticeable also in the metallic 
phase |jo| B3, |95| . The second pool of states determines 
the conductivity, so that properties of these states are 
described by the scaling theory of MIT. Accordingly, the 
corresponding localisation radius £ is rather controlled 
by interference of multi-scattering processes than by the 
attractive potential of a single impurity. Thus, £ of these 
weakly localised states is significantly larger than as, and 
diverges on approaching the MIT from the insulator side. 
It is worth noting that such a two-fluid model is consis- 
tent with a.c. conductivity studies Q, which show the 
coexistence of weakly and strongly localised states near 
the MIT in (Ga,Mn)As. Furthermore, the merging of im- 
purity and band states in this range is substantiated by 
angle-resolved photoemission spectra in the same system 




In order to tell the dominant mechanism accounting 
for the existence of long-range spin order in ferromag- 
netic semiconductors it is instructive to trace the evolu- 
tion of their magnetic properties on crossing the MIT. 
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Remarkably, in contrast to rather strong changes of re- 
sistivity, the evolution of magnetic properties is gradual. 
This substantiates the notion that thermodynamic prop- 
erties do not exhibit any critical behaviour at MIT as 
they are insensitive to large-scale characteristics of the 
wave functions. Importantly, the values of Curie tem- 
perature are found to grow with the degree of the ma- 
terial metallicity jl(| ^8|. Moreover, the exam- 
ination of magnetization as a function of temperature 
and magnetic field indicates that virtually all Mn spins 
contribute to ferromagnetic order in the most metallic 
samples |§ ||, M ^ However, on crossing the 
MIT (by lowering x) , the relative concentration of ferro- 
magnetically coupled spins decreases substantially. Ac- 
cording to XMCD results 0, about 10% of Mn spins 
is involved in ferromagnetism of Gai-^Mn^As with x = 
2%. Also ferromagnetic resonance studies Q and di- 
rect magnetization measurements demonstrate that only 
a part of spins contribute to spontaneous magnetization, 
while the alignment process of the remaining moments 
occurs according to a Brillouin function for a weakly in- 
teracting spin system J9l[ . Remarkably, the anomalous 
Hall effect reveals clearly the presence of the first compo- 
nent but hardly points to the existence of any loose spins 

@ n. 

The above findings indicate that Mn spins in the re- 
gions visited by itinerant holes are coupled ferromagnet- 
ically. These holes set long-range ferromagnetic correla- 
tion between Mn spins, including those contributing to 
BMP that are formed around singly occupied local states. 
Obviously, the ferromagnetic portion of the material, and 
thus the magnitude of spontaneous magnetization, grows 
with the dopant concentration, attaining 100% in the 
metallic phase. Such a trend is confirmed by the avail- 
able data, as discussed above. Thus, the delocalised or 
weakly localised holes are responsible for ferromagnetic 
correlation in both (Ga,Mn)As and (Zn,Mn)Te:N [p7|. 
At the same time, mechanisms that involve strongly lo- 
calised states, such as excitations from impurity levels or 
a direct coupling between BMP, appear to be of lesser 
importance. 

According to the two-fluid scenario referred to above, 
the BMP are present on the both sides of the metal- 
insulator transition (MIT) |95|j . As already mentioned, 
the coupling between the BMP appears to be ferromag- 
netic JS(| . Since Tq is proportional to the square of spin 
vector length, the weight of the BMP contribution may 
exceed their relative concentration. To gain the Coulomb 
energy, the BMP are preferentially formed around close 
pairs of ionised acceptors. In the case of III-V mate- 
rials, one hole localised at two Mn ions generates, via 



Zener's double exchange 70 , a strong ferromagnetic cou- 
pling that overcompensates the intrinsic superexchangc 
antifcrromagnetic interaction. In contrast, in II- VI com- 
pounds, for which the acceptor cores do not carry any 
spin and the degree of compensation by donors is low, 



BMP are not preferentially form around the Mn pairs. 
Accordingly, the pairs of the close Mn spins remain anti- 
ferromagnetically aligned, even in p-type samples. This 
diminishes the effective concentration of the Mn spins 
contributing to the ferromagnetic order, x e g < x, and 
decreases the effective Curie- Weiss temperature by T^p, 
both parameters known from magnetization studies on 
undopcd II- VI DMS. The weakening of the ferromag- 
netism by the superexchange in the case of II- VI com- 
pounds, and the virtual absence of the corresponding ef- 
fect in III-V materials, where ?Ap ~ and x e ^ w x, con- 
stitutes the important difference between these two fam- 
ilies of magnetic semiconductors. This fact is taken into 
account when evaluating the Mn contribution to the total 
free energy of particular systems. In worth adding that 
recent Monte Carlo simulations j9|] carried out start- 
ing from an impurity band model, provides an additional 
support for the two fluid scenario. 



Mean-field Zener model and its limitations 

In terms of the mean-field Zener model of the carrier 
mediated ferromagnetic interactions pq, 



10C, 101 



the equilibrium magnetization, and thus the Tq is de- 
termined by the minimum of the Ginzburg-Landau free- 
energy functional F[M(r)] of the system, where M(r) 
is the local magnetization of the localised spins. This 
is a rather versatile scheme, to which car rier correlation 
and confinement fl5| , [Tsj , 39|, 101 , 102 \, k ■ p and spin- 
orbit couplings |l6|, |27|, p8[ 103 1 as well as weak disor- 



der and antiferromagnetic interactions [lj, |26| can 



be introduced in a controlled way, and within which the 
quantitative comparison of experimental and theoretical 
results is possible Jl(| |2(], |2^, |l 0-4f| . In its general formu- 
lation, the model allows for non-uniform ground states 
(M(q > 0) ^ 0), such spin-density waves or non-collinear 
(counted) magnetic structures 

In theory developed by the present authors and co- 
workers [^6[ 27, [28fl , the hole contribution to F is com- 
puted by diagonalizing the 6x6 Kohn-Luttinger k ■ p- 
matrix containing the p-d exchange contribution, and by 
the subsequent computation of the partition function Z, 
F c = —ksTlnZ. In the case of a strongly degenerate 
Fermi liquid, — £/p|/fcgT 3> 1, F c can be replaced by the 
ground state energy. However, this approximation would 
severely overestimate T(j in materials with low hole den- 
sities or with high Tq values and, therefore, has not been 
employed for such situations |2C], ^?], ^8). The model is 
developed for both zinc-blende and wurzite semiconduc- 
tors, takes the effects of the spin-orbit interaction into 
account, and allows for the presence of both biaxial strain 
and quantizing magnetic field. The enhancement of the 
tendency towards ferromagnetism by the carrier-carrier 
exchange interactions is described by the Fermi liquid pa- 
rameter Ap. The value A-p = 1.2 was evaluated within 
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the LSDA for the 3D carrier liquid of the relevant density 
[ 101 1 . Importantly, the formalism has been extended by 



Konig et al. [105, 106 1 , who evaluated the dependence 
of the carrier free energy on the wave vector q. This 
dependence provides information on magnetic stiffness, 
which together with magnetic anisotropy, determine the 



netic domains 104 



dispersion of spin waves [106] and the structure of mag 



It is, of course, important to comment explicitly main 
approximations behind the Zener model. Obviously, this 
model may not be applicable to DMS containing mag- 
netic impurities other than Mn, in which d levels reside 
in the band gap and correlation energy U is relatively 
small. In the Mn-based DMS, however, the magnetic 
constituent appears to be well described in terms of the 
the Anderson impurity model and, thus, such DMS can 
be regarded as charge transfer insulators. This means 
that the spins are localised, so that the d electrons do 
not participate in charge transport. There exist, how- 
ever, quantum fluctuations (hybridization) between the 
p and d orbitals, which result in the Kondo-likc inter- 
action that couples the spin and the carrier parts of the 
free energy. It should be emphasized that results of the 
LSDA computations |7j| point to a rather large weight of 
the d electrons at the Fermi level and, hence, imply the 
ferromagnetism to be driven by the mechanism of double 
exchange. This controversy is perhaps the most intrigu- 
ing open issue in theory of magnetic semiconductors. 

The use of the spin, not carrier magnetization, as the 
order parameter in the Zener model is a consequence of 
the adiabatic approximation: spin dynamics is assumed 
to be much slower than that of the carriers. This approx- 
imation may break down if, for instance, a characteristic 
energy of spin-spin antiferromagnetic interactions would 
become greater than Tq. Furthermore, the order param- 
eter can be regarded as continuous as long as the concen- 
tration of the spins is much larger than that of the holes, 
x e gNo 3> p. Importantly, in this range, the mean-field 
value of the ordering temperature T(j(q) deduced from 
the Zener and the RKKY model are identical, indepen- 
dently of microscopic spin distribution [ p6[ . 

However, in the opposite limit, x e gNo < p, impor- 
tant changes in the hole response function occur at the 
length scale of a mean distance between the localised 
spins. Accordingly, the description of spin magnetization 
by the continuous-medium approximation, which consti- 
tutes the basis of the Zener model, ceases to be valid. 
In contrast, the RKKY model is a good starting point 
in this regime Jl6| , as it provides the dependence of the 
interaction energy of particular spin pairs as a function 
of their distance. This makes it possible to evaluate the 
system energy for a given distribution of the localised 
spins. The resulting competition between the ferromag- 
netic and antiferromagnetic interactions is expected to 
grow with p/x e ffNQ, and may lead, through non-collinear 
(counted) spin arrangement [107, 108], to the spin-glass 



phase [ 109 1 . Alternatively, with decreasing x at given p, 
the Kondo effect may show up |26| . 



It is interesting to note that the role of thermodynamic 
fluctuations of magnetization, that is the inaccuracy of 
the mean-field approximation (MFA), grows also with 
p/x e ffNo. It is well known that the MFA results are 
exact, also in low-dimensional systems, if the range of 



ferromagnetic interactions is infinite [ 1 1 1 . The decay of 
the strength of the carrier-mediated exchange interaction 
with the distance r between two Mn spins is described by 
the RKKY oscillatory function. At small r, the interac- 
tion is ferromagnetic, and then changes sign at r = 1.2r c , 
where r c is an average distance between the carriers that 
mediate the spin-spin coupling. This implies that the 
MFA is valid quantitatively at p << x e qNo, a conclu- 
sion consistent with the estimate of Tq taking the spin 
wave excitations into account [105]. Actually, however, 
the range of validity of the MFA is s i gnifi cantly larger 
[ 107 1 than that initially proposed [105, 111], as the mag- 
nitudes of spin stiffness evaluated within the 6x6 Lut- 
tinger model is much greater [106] than those obtained 
for a simple parabolic band 1 105 . The dynamic mean- 



field approach [112| and, especially, hybrid Monte-Carlo 
algorithms 113 have potential to sheet some light on 
ferromagnetism in the regime beyond the validity of the 
MFA. 



Finally, we address the question about the role played 
by disorder. Since thermodynamic properties of the car- 
rier liquid are relatively weakly perturbed by scattering, 
its effect on the hole-mediated ferromagnetism can be 
neglected to a first approximation. When disorder grows 
and the MIT is approached, the mean free path becomes 
comparable to the inverse Fermi wave vector. Within the 
Zener model, the effect of the finite mean free path l e can 
be described by scattering broadening of the density of 
states [|~8[ [26| , which reduces Tq. Technically, this fol- 
lows from the averaging of the free energy over possible 
impurity distributions. Equivalently, diffusive character 
of carrier transport leads to exponential dumping of the 
RKKY interactions at distances longer than l e . How- 
ever, large fluctuations in the carrier distribution have to 
be taken into account at criticality and on the insulator 
side of the MIT. As already argued, the magnetization 
is small in areas which are not visited by the carriers, 
whereas its magnitude is large there, where delocalised 
or weakly localised carriers reside. Remarkably, this en- 
hances TV; over the value expected for the average carrier 
density Q but reduces the magnitude of sample-average 
spontaneous magnetization (If], |5 §§. The interplay 
between Anderson-Mott localisation, Stoner magnetism, 
and carrier-mediated spin-spin interaction is certainly an 
appealing area for future research. Recent theoretical 



works in this direction are indeed encouraging 1 1 14 ] 
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COMPARISON OF THE ZENER MODEL TO 
SELECTED EXPERIMENTAL RESULTS 

Magnetic circular dichroism in (Ga,Mn)As 

Within the Zener model, the strength of the ferromag- 
netic spin-spin interaction is controlled by the k-p param- 
eters of the host semiconductor and by the magnitude of 
the spin-dependent coupling between the effective mass 
carriers and localised spins. In the case of II- VI DMS, de- 
tailed information on the exchange- induced spin-splitting 
of the bands, and thus on the coupling between the ef- 
fective mass electrons and the localised spins has been 
obtained from magnetooptical studies 1 52 , |5j| . A similar 
work on (Ga,Mn)As [j]]], |l| led to a number of 

surprises. The most striking was the opposite order of 
the absorption edges corresponding to the two circular 
photon polarizations in (Ga,Mn)As comparing to II- VI 
materials. This behaviour of circular magnetic dichroism 
(MCD) suggested the opposite order of the exchange-split 
spin subbands, and thus a different origin of the sp-d in- 
teraction in these two families of DMS. A new light on the 
issue was shed by studies of photoluminescence (PL) and 
its excitation spectra (PLE) in p-type (Cd,Mn)Te quan- 
tum wells pl|. As shown schematically in Fig. 3, the 
reversal of the order of PLE edges corresponding to the 
two circular polarizations results from the Moss-Burstein 
effect, that is from the shifts of the absorption edges as- 
sociated with the empty portion of the valence subbands 
in the p-type material. This model was subsequently 
applied to interpret qualitative ly th e magnetoabsorption 
data for metallic (Ga,Mn)As |117|] . More recently, the 
theory was extended by taking into account the effect of 
scattering-induced mixing of k states |58). As shown in 
Fig. 4, this approach explains the slop of the absorption 
edge as well as its field-induced splitting assuming the 
value of the p-d exchange energy (3 No = — 1 cV. 

Surprisingly, however, the anomalous sign of the MCD 
was present also in non-metallic (Ga,Mn)As, in which 
EPR signal from occupied Mn acceptors was seen |36|| . 
It has, therefore, been suggested that the exchange inter- 
action between photo- and bound-holes is responsible for 



the anomalous sign of the MCD in those cases [117]. The 
presence of such a strong exchange mechanism is rather 
puzzling, and it should be seen in non-magnetic p-type 
semiconductors. At the same time, according to our two- 
fluid model, the co-existence of strongly and weakly lo- 
calised holes is actually expected on the both sides of the 
MIT. Since the Moss-Burstein effect operates for inter- 
band optical transitions involving weakly localised states, 
it leads to the sign reversal of the MCD, also on the in- 
sulating side of the MIT. 

Another striking property of the MCD is a different 
temperature dependence of the normalized MCD at low 
and high photon energies in ferromagnetic (Ga,Mn)As 
[ L LS | . This observation was taken as an evidence for the 



m 
i— i 



OOe 



50 Oe 



150 Oe 



PLE 



j — i 




50 kOe 



1640 1660 1680 1700 
Energy (meV) 



FIG. 3: Photoluminescence excitation spectra (PLE), that is 
the photoluminescence (PL) intensity as a function of the ex- 
citation photon energy intensity, for a + (solid lines) and a~ 
(dotted lines) circular polarizations at selected values of the 
magnetic field in a modulation-doped p-type quantum well 
of Cdo.976Mno.o24Te at 2 K. The photoluminescence was col- 
lected in a + polarization at energies marked by the narrowest 
features. The sharp maximum (vertical arrow) and step-like 
form (horizontal arrow) correspond to quasi-free exciton and 
transitions starting at the Fermi level, respectively. Note re- 
verse ordering of transition energies at a + and a~ for PL 
and PLE (the latter is equivalent to optical absorption). The 
band arrangement at 150 Oe is sketched in the inset (after 
Ref. Pi). 



presence of two spectrally distinct contributions to op- 
tical absorption [118|. A quantitative computation of 
MCD spectra has recently been undertaken pi|. The 
theoretical results demonstrate that because of the Moss- 
Burstein effect, the magnetization-induced splitting of 
the bands leads to a large energy difference between the 
positions of the absorption edges corresponding to the 
two opposite circular polarizations. This causes an un- 
usual dependence of the low-energy onset of MCD on 
magnetization, and thus on temperature. These consid- 
erations lead to a quantitative agreement with the ex- 
perimental findings, provided that the actual hole disper- 
sion and wave functions are taken for the computation of 
MCD. 
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FIG. 4: Transmission of Gao.96sMno.rj32As film for two cir- 
cular light polarizations in the Faraday configuration in the 
absence of the magnetic field (data shifted up for clarity) and 
in 5 T at 2 K (points) 117 1. Solid lines are calculated for 



the hole concentration p — 7 x 10 19 cm -3 , exchange energy 
No/3 = —1 eV, and allowing for scattering-induced breaking 
of the k selection rules [681 . 



In conclusion, a.c. conductivity in the far infrared 
as well as to photoemission |59| |6(| and XMCD 
|6^ | in the range of high photon energies, together 
with magnetooptical characteristics discussed here, are 
consistent with the picture of electronic states advocated 
for (Ga,Mn)As in the previous section. In particular, 
(Ga,Mn)As exhibits properties generic to doped semi- 
conductors in the vicinity of the metal-insulator tran- 
sition. Furthermore, no experimental results have so 
far been collected in favor of the valence band splitting 
corresponding to the p-d exchange integral as high as 
/3Nq —4 eV, as suggested by ab initio computations 
within the LSDA Ifl. 



Curie temperature and spontaneous magnetization 
in p-type and n-type Mn-based DMS 

In order to compare theoretical expectations con- 
cerning Tq to experimental results, it is convenient 
to introduce the normalized ferromagnetic temperature 
Tp/x e g = (Tq +2Ap)/a; e ff which, within the mean-field 
Zener model, should not depend on the Mn concentration 
x. We recall that x e q < x and T^p > take into account 
the presence of antiferromagnetic exchange interaction in 
II- VI DMS. Figure 5 presents Tp /z» ff for Gai-^Mn^As 



16|, 120 1, and quantum 
as a function the Fermi 



@ H [m§, p-Zm^Mn,!^ 
wells of p-Cdi-zMn^Te 
wavevector determined from the value of the hole concen- 
tration p assuming the Fermi sphere to be isotropic. The 
hole concentration was deduced either from the Hall re- 
sistance jl6, 119 or form the Moss-Burstein shift |l^, p0| . 



FIG. 5: Experimental (symbols) and calculated (lines) nor- 
malized ferromagnetic temperature, Tp/10 2 :r e g-, versus the 
wav e ve ctor at the Fermi level for Ga j -^M n^As (triangle) 



c Mn x Te:N (squares) [JL6 



12J, Zni-zMn^TeiP 
Mn^Te (circles) 



(sta r'l JL2C| ], and quantum well of p-Cdi- 
|[18| , p3| . Solid lines: Zener and 6x6 Luttinger model for 
the 3D Q) and 2D case dotted line: the RKKY and 
6x6 Luttinger model for x e g = 0.015, taking into account 
the effect of the antiferromagnetic interactions on statistical 
distribution of unpaired Mn spins [ |l6[ . 



We have to emphasize that because of a contribution 
from the anomalous (spin) Hall effect and a non-uniform 
hole distribution in thin layers, the evaluation of the hole 
concentration is by no means straightforward in magnetic 
semiconductors . 

A number of important conclusions emerges from 
the comparison of experimental and theoretical results. 
First, the theory ^ || with (3N = -1.2 eV p| ex- 
plains the large magnitude of Tq = 110 K jl4], ]8^] for 
Gao. 947Mno.Q53As containing 3.5 x 10 20 holes per cm 3 
[mj. Second, in contrast to Mn-based III-V DMS, it is 
essential — in the case of II- VI materials — to take into ac- 



count the presence of magnetically inert nearest-neighbor 
Mn pairs. Such pairs not only lowers the effective Mn 
concentration but also make the antiferromagnetic por- 
tion of the RKKY interaction to become more significant, 
which lowers Tq at large p (dotted line) p6j . Third, 
scaling theory of electronic states near the MIT, dis- 
cussed in the previous section, makes it possible to ex- 
plain the presence of ferromagnetism on the both sides 
of the MIT, and a non-critical evolution of Tq across the 
critical point, a behaviour observed in both (Ga,Mn)As 
p3 and p-(Zn,Mn)Te p|, pM Imf. Importantly, in 



agreement with this scenario, a ferromagnetic component 
of the material increases with the hole concentration [jl6| . 
However, as already explained, because of RKKY oscil- 
lations, the present theory (solid lines) overestimates the 
magnitude of Tq if the hole and Mn concentrations be- 
comes comparable, p ~ x e ff-/Vrj. In contrast, the actual 
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value of Tq can be larger than that of Fig. 5 on the insu- 
lator side of the MIT because of the fluctuations in the 
carrier density and a non-linear contribution of spin clus- 
ters within the BMP. This might be the reason for a rel- 
atively high values of Tq detected recently in (In,Mn)As 
[ 122 1 . Last but not least, according to Fig. 5, in mod- 
ulation doped (Cd,Mn)Te/(Cd,Mg,Zn)Te:N heterostruc- 
tures, due to the enhanced density-of-states (DOS) at low 
energies in 2D systems and reduced localisation by the 
distant nitrogen acceptors, the ferromagnetic transition 
occurs even for the densities of the hole liquid as low as 
10 17 cm~ 3 [[l5| [ll|. The present theory describes cor- 
rectly the dependence of Tq on x and p, provided that 
the exchange interaction between the holes and disorder 
broadening of DOS are taken into account [[l8j [20[ . 

Two effects appear to account for the greater Tq values 
in p-(Ga,Mn)As than in p-(Zn,Mn)Te at given p and x. 
First is the smaller magnitude of the spin-orbit splitting 
between the Tg and T7 bands in arsenides, A D = 0.34 eV, 
in comparison to that of tellurides, A D = 0.91 eV. 
Once the Fermi energy E-p approaches the band, the 
density-of-states effective mass increases, and the reduc- 
tion of the carrier spin susceptibility by the spin-orbit 
interaction is diminished. The computed value of Tq 
for p = 3 x 10 20 cm~ 3 is greater by a factor of four in 
(Ga,Mn)As than that evaluated in the limit A Q >> E-p. 
The other difference between the two materials is the de- 
structive effect of antifcrromagnetic interactions, which 
operate in II- VI compounds but are of minor importance 
in III-V materials, as argued in the previous section. 

Since in semiconductors the magnitude of exchange 
splitting of bands is comparable to the Fermi energy, 
the growth of spontaneous magnetization M s on lowering 
temperature deviates from the Brillouin-typc behaviour 
even in the mean-filed approximation [^5[ [l(| [H| |s) . 
Because of extremely low hole densities, the effect is par- 
ticularly well visible in the case of modulation-doped 
(Cd,Mn)Te/(Cd,Mg,Zn)Te:N heterostructures || [18). 
Furthermore, the contribution from carrier magnetic mo- 
ments and disorder-induced fluctuations in the carrier 
distribution constitute other sources of corrections to 
the Brillouin-type dependence. However, with the in- 
crease of the hole concentration, M S (T) evaluated within 
the mean-field Zener model tends toward the Brillouin 
function for material parameters of (Ga,Mn)As pgfl , an 
expectation corroborated by the experimental results 

Finally, we turn to materials, in which the Fermi level 
resides in the s-type conduction band. Because of small 
s-d exchange energy and low density of states, no ferro- 
magnetism is expected above 1 K in such a case pa . Ex- 



Effects of strain 

Both hydrostatic and axial strain affect the valence 
band, and thus alter the magnitude of the density of 
states and Tq. Quantitatively, however, the effect is 
evaluated to be small pij] . There exists another mecha- 
nism by which strain may affect Tq. It is presently well 
known that the upper limit of the achievable carrier con- 
centration is controlled by pinning of the Fermi level by 
impurity or defect states in virtually all compound semi- 
conductors. Since the energies of such states in respect 
to bands vary strongly with the bond length, the hole 
concentration and thus Tq will depend on strain. 

Apart from Tq and spontaneous magnetization M s , 
it is interesting to consider means making it possible 
to tailor magnetic anisotropy, and thus the direction of 
the spontaneous magnetization, the coercive force, the 
switching field, the domain structure. Already early 
studies of the ferro magnetic phase in Ini-^Mn^As | 12C | 



and Gai-^Mn^As [127 



128 



pcrimental results for n-(In,Mn)As [ 123 1 , (Ga,Mn)As:Sn 
[ 124] , and (Zn,Mn)0:Al ^2(J confirm this expectation. 
With this in mind, the presence of indications of ferro- 



magnetism in n-type (Ga,Mn)N [ 125 1 is challenging 



demonstrated the existence 
of sizable magnetic anisotropy. Magnetic anisotropy is 
usually associated with the interaction between spin and 
orbital degrees of freedom of the d-electrons. Accord- 
ing to the model advocated here, these electrons are in 
the d 5 configuration. For such a case the orbital mo- 
mentum L = 0, so that effects stemming from the spin- 
orbit coupling are expected to be rather weak. It has, 
however, been noted that the interaction between the 
localised spins is mediated by the holes that have a non- 
zero orbital momentum p7fl . An important aspect of 
the Zener model is that it does take into account the 
anisotropy of the carrier-mediated exchange interaction 
associated with the spin-orbit coupling in the host mate- 
103 1, an effect difficult to include within the 



rial {27J |28| 

standard approach to the RKKY interaction 10? ]. 

A detail numerical analysis of anisotropy energies has 
been carried out for a num ber of experimentally im- 
portant cases |2?j, |28|, 103 . In particular, the cubic 
anisotropy as well as uniaxial anisotropy under biaxial 
strain have been studied as a function of the hole con- 
centration p. The computation indicates that for the 
parameters of Gai-^Mn^As films grown along the [001] 
direction, the spontaneous magnetization M lies in the 
(001) plane, and the easy axis is directed along [100] or 
along [110] (or equivalent) crystal axis depending on the 
degree of the occupation the hole subbands as well as on 
their mixing by the p-d and k ■ p interactions. As a re- 
sult, the easy axis fluctuates between [100] and [110] as a 
function of p, the preferred direction for typical hole con- 
centrations being [110]. The magnitude of the external 
magnetic field H cu that aligns M along the hard direc- 
tion in the (001) plane is evaluated to be up to 0.2 T p8| . 
Since, however, the orientation of the easy axis changes 
rapidly with p and M, disorder-which leads to broaden- 
ing of hole subbands-will presumably diminish the actual 
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magnitude of magnetic anisotropy. The field ^ H CU de- 
termines also the magnitude of the switching field, which 
could be observed in microstructures containing only a 
single domain. In macroscopic films, however, smaller 
values of the coercive field fi H c are expected as actually 
observed: typically fi H c — 4 mT for the magnetic field 
along the easy axis in the (001) plane in Gai_ K Mn x As 
[12§. 

It can be expected that strain engineering can ef- 
ficiently control magnetic properties resulting from 
the hole-mediated exchange. Indeed, sizable lattice- 
mismatch driven by biaxial strain is known to exist 
in semiconductor layers. In some cases, particularly if 
epitaxy occurs at appropriately low temperatures, such 
strain can persist even beyond the critical thickness due 
to relatively high barriers for the formation of misfit dis- 
locations. It has been found that the biaxial strain leads 
to uniaxial anisotropy, whose magnitude can be much 
greater than that resulting from either cubic anisotropy 
or stray fields. As shown in Fig. 6, for the experimen- 
tally relevant values of p and M, the easy axis is pre- 
dicted to be oriented along [001] direction for the ten- 
sile strain, whereas it should reside in the (001) plane 
for the c ase o f unstrained or compressively strained films 
102 L_ This is corroborated by the experimen- 



tal study [127, 128 1, in which either (In,Ga)As or GaAs 
substrate was employed to impose tensile or compres- 
sive strain in Gai-^Mn^As, respectively. In particu- 
lar, for the Gao.965Mno.o35As film on GaAs, for which 
(-xx = -0.24%, the anisotropy field fJ, H un = 0.4 ± 0.1 T 
is observed [127 , 



in quantitative agreement with the 
theoretical results of Fig. 6. This field is about two or- 
ders of magnitude greater than that evaluated from the 
extrapolation of ESR data on single- ion anisotropy at low 
x [ p_29| , a result confirming the dominant contribution of 
the holes to the magnitude of H un . Though no theoreti- 
cal computations have been performed for Ini-^Mn^As, 
the qualitatively similar effect of biaxial strain is ex- 
pected, in agreement with the early experimental results 



[126| 



It worth noting that similarly to strain, also confine- 
ment of the holes affects the magnetic anisotropy - in 
accord with the theoretical model, the easy axis is ori- 
ented along the growth direction in the ferromagnetic 
p-(Cd,Mn)Te quantum wells g @. 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

Biaxial strain e (%) 



FIG. 6: Computed minimum magnetic field H un necessary 
to align the saturation value of magnetization M 3 along the 
hard axis as a function of biaxial strain component e xx for 
two values of the hole concentrations in Gao.95Mno.05 As [Q. 
The symbol [100] — > [001] means that the easy axis is along 
[100], so that H un is applied along [001]. 
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FIG. 7: Temperature dependence of the width of domain 
stripes as measured for the Ga0.957Mn0.043As fi lm w ith the 
easy axis along the growth direction (full squares) |l3Cj | . Com- 
puted domain width is shown by the solid line. The dashed 
fine is computed assuming that the parameter A c (the ratio 
of the dom ain w all and stray field energies) is by a factor of 
1.8 greater fo3]. 



Domain structure 

Recently, the structure of magnetic domains in 
Gai-ajMn^As under tensile strain has been determined 
by micro-Hall probe imaging ]l30fl . The regions with 
magnetization oriented along the [001] and [00 ] easy 
axis form alternating stripes extending in the [110] di- 
rection. This indicates, for either Bloch or Neel domain 
walls, that the in-plane easy axis is rather along [110] 



than along [100] directions, a conclusion consistent with 
the theoretical expectation for in-plane (cubic) magnetic 
anisotropy presented above. As shown in Fig. 7, the ex- 
perimentally determined stripe width is W = 1.5 /xm at 
5 K for 0.2 fi film of Gao.957Mno.043 As on Gao.s4Ino.i6As, 
for which tensile strain of e xx — 0.9% is expected. 



According to micromagnetic theory [ 131 j , W is de- 
termined by the dimensionless parameter A c , which is 
given by the ratio of the domain wall and stray field 
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energies. The former is proportional to M 2 , whereas 
the latter scales with the product of anisotropy energy 
K u [^8|, 103 1 and magnetic stiffness A ]106[ . Figure 7 
presents the calculated values of W(T) Jl04[ in compar- 



ison to the experimental data for (Ga,Mn)As [130]. The 
material parameters collected in Ref. |28| , and employed 
to generate theoretical results of Figs. 5 and 6, have been 
adopted. Furthermore, in order to establish the sensi- 
tivity of the theoretical results to the parameter values, 
the results calculated for a value of A c 1.8 times larger 
are included as well. The computed value for low tem- 
peratures, W = 1.1 /im, compares favorably with the 
experimental finding, W — 1.5 /im. However, the model 
predicts much weaker temperature dependence of the do- 
main w idth W than observed experimentally, which are 
linked [ 1 04 1 to critical fluctuations, disregarded in the 
mean-field approach. 



TOWARDS HIGH-TEMPERATURE 
FERROMAGNETIC SEMICONDUCTORS 

Mn-based tetrahedrally coordinated DMS 

In view of the general agreement between experiment 
and theory for Tq and the magnetic anisotropy, it is 
tempting to extend the model for material systems that 
might be suitable for fabrication of functional ferro- 
magnetic semiconductors. For instance, the model sug- 
gests immediately that Tq values above 300 K could 
be achieved in Gao.gMn .iAs, if such a large value of 
x would be accompanied by a corresponding increase in 
the hole concentration. Figure 8 presents the values of Tq 
computed for various tetrahedrally coordinated semicon- 
ductors containing 5% of Mn per cation and 3.5 x 10 20 
holes per cm 3 |27|]. In addition to adopting the tabu- 
lated values of the band structure parameters, the same 
value of [3 = [3{G&i- x Mn x ks] for all group IV and III-V 
compounds was assumed, which results in an increase of 
\PNq\ ~ a~ 3 , where a is the lattice constant, a trend 
known to be obeyed within the II- VI family of magnetic 
semiconductors. For the employed parameters, the mag- 
nitude of Tq for the cubic GaN is by 6% greater than 
that computed for the wurzite structure. 

The data (Fig. 8) demonstrate that there is much room 
for a further increase of Tq in p-type magnetic semicon- 
ductors. In particular, a general tendency for greater 
Tq values in the case of lighter elements stems from the 
corresponding increase in p-d hybridization and the re- 
duction of the spin-orbit coupling. It can be expected 
that this tendency is not altered by the uncertainties in 
the values of the relevant parameters. Indeed, the re- 
sults of Fig. 8 have triggered a considerable fabrication 
effort, which is bringing a number of striking develop- 
ments. In particular Gei^^Mn^ was found to be p-type 
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FIG. 8: Computed values of the Curie temperature Tq for 
various p-type semiconductors containing 5% of Mn per cation 
f2.5% per atom) and 3.5 x 10 20 holes per cm 3 (after Refs. 



higher magnitude of Tq, in excess of 200 K, was found in 
(Ga,Mn)P:C |133f| , a value consistent with the prediction 
of Fi g. 8 |pJJ p8| for the employed Mn content x — 9.4% 
[ 133 1 . Moreover, there are indications of ferromagnetism 
in (Ga,Mn)N Jl34|, HI [l3^| in some cases with 



T c near or above 300 K [|135j [136[ |137| . While further 



and to exhibit Tq in the excess of 100 K [132]. An even 



experimental work is needed to elucidate a possible role 
of precipitates of various Mn-Ga and Mn-N compounds, 
the highest value s ugg ested up to now, Tq = 940 K 
for Gao.91Mno.09N [ 136 1 , is consistent with the expecta- 
tions of the mean-field Zener model ]^9| . In more general 
terms, within this model, large magnitudes of Tq result 
from the combine effect of the large on-site p-d exchange 
interaction and the efficient transfer of spin information 
owing to a relatively large extend of the p-wave functions. 

Another way pursed recently consists of searching for 
the high temperature ferromagnetism in compounds con- 
taining combination of elements from both II- VI and III- 
V materials. Layered Gai-^Mn^Se [ 138 1 and particu- 
larly chalcopyrite Cdi__ a; Mn 2 ;GeP2 ]139[ are starting to 
provide encouraging results. Since, in general, III-V com- 
pounds can easier be doped by impurities that are elec- 
trically active, whereas II- VI materials support greater 
concentration of transition metals, a suggestion has been 
put forward to g row magnetic III-V/II-VI short period 
superlattice |140| , in which a charge transfer to the mag- 
netic layers will increase Tq. 



Beyond Mn-based compounds 

The highest spin and the associated large magnitude 
of the on-site correlation energy U account for the di- 
valent character of the Mn atoms in a large variety of 
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environments. This results, in particular, in a large sol- 
ubility of Mn in II- VI materials and its acceptor charac- 
ter in III-V compounds. A question arises about ferro- 
magnetic properties of semiconductors containing other 
magnetic components. Quite generally, carrier-mediated 
exchange mechanism is efficient if both concentration of 
holes and p-d exchange energy ft No are appropriately 
large. According to Fig. 1, for many transition metals 
and hosts, the d states lie in the gap. Hence, to trig- 
ger carrier-mediated exchange interactions, co-doping by 
shallow impurities is necessary. However, shallow accep- 
tors will produce band holes only if the d N+1 /d N donor 
level resides at lower energies, that is in the valence band. 
Furthermore, except for some resonant situations, the 
magnitudes of \(3Nq\ are small for N < 5 because of oppo- 
site signs of two relevant contributions [141|. With these 
arguments in mind, (Ga,Fe)N and (Ga,Co)N co-doped 
with acceptor impurities appear as promising materials 

urn- 

There are, however, other possible roads to high tem- 
perature ferromagnetism. Since for configurations other 
than d 5 , the magnitude of U is relatively small, the 
Hubbard-Mott transition, and the associated ferromag- 
netism may appear at increasing density of the magnetic 
constituent. To our knowledge no such effects have so 
far been detected in materials in question. Furthermore, 
for some magnetic ions and hosts, either d N /d N+1 or 
d N /d N ~ x level lies below the expected energy of shal- 
low donor or above energy of shallow acceptor, respec- 
tively. In such a situation, the co-doping by shallow im- 
purities will rather affect the occupancy of the d levels 
than introduce carriers into the bands, a known semi- 
insulating character of semiconductors containing spe- 
cific transition metals. At appropriately large concen- 
trations of both magnetic ions and additional impurities, 
the carriers trapped on d levels may start to transfer, via 
Zener's double exchange mechanism |f7(i|| , magnetic infor- 
mation between localised spins. Indications of ferromag- 
netism in the vicinity of room temperature discovered in 
n-(Zn,Co)0 |U§ and n-(Zn,V)0 (m| may represent, ac- 
cording to Fig. 1, such a case. Actually, a simultaneous 
doping by two different magnetic impurities may serve 
to control both spin and charge states as well as the 
carrier concentration in the specific bands. The room- 
temperature ferromagnetism detected after doping a few 
percent of Co to nonmagnetic TiC>2 [ 144 1 may belong to 
this category. In particular, the hybridized d levels of Ti 
and Co may constitute the efficient channel for transmis- 
sion of magnetic information between the Co spins. 

Finally, one should recall the existence of, e.g., ferro- 
magnetic europium chalcogcnides and chromium spinels. 
In those compounds, ferromagnetism is not mediated by 
carrier transport. With no doubt, the availability of 
intrinsic and n-type tetrahedrally-coordinated ferromag- 
netic compounds would enlarge considerably the impact 
of semiconductor electronics. Actually, a theoretical sug- 



gestion has been made 1 74, 75 that superexchange in 
Cr-based and V-based II- VI compounds can lead to a 
ferromagnetic order. Desired material properties, such as 
divergent magnetic susceptibility and spontaneous mag- 
netization, can also be achieved in the case of a strong 
antifcrromagnetic super-exchange interaction. The idea 



here [ 145 1 is to synthesize a ferrimagnetic system that 
would consist of antiferromagnetically coupled alternat- 
ing layers containing different magnetic cations, e.g., Mn 
and Co. 

It is also tempting to consider the performance of DMS 
containing ions from other groups of magnetic elements. 
Current works on Sii-zCe^ [ 146 1 and on the insulator-to- 
metal transition in amorphous Sii-^Gda, [147 1 alloys can 



be quoted in this context. However, because of a weak hy- 
bridization between 4f and band states (exploited in Er- 
doped emitters) , no high-temperature ferromagnetism is 
expected in rare-earth-based systems. More promising in 
this respect are materials containing 4d or 5f ions — there 
exists already a preliminary report on spin-spin interac- 
tions in undoped Pbi-zU^Te [148]. Actually, in view of 
indications of room-temperature weak ferromagnetism in 
(Ca,La)B6 Jl49| and polimerized rh-Cgo [ 150 , which are 
built up of non-magnetic constituents, the spectrum of 
possibilities appears as unlimited. 

The above list of achievements and prospects clearly 
shows that searches for high temperature ferromagnetic 
semiconductors have evolved into a broad field of materi- 
als science. In addition to work on design and synthesis 
of new systems, a considerable effort will certainly be 
devoted to control contributions from ferromagnetic or 
ferrimagnetic precipitates and inclusions in the materi- 
als available already now. On theoretical side, the in- 
terplay between Anderson-Mott localisation, disordered 
Stoner magnetism, and carrier-mediated ferromagnetism 
will attract a considerable attention. With no doubt we 
will witness many unforeseen developments in the field 
of ferromagnetic semiconductors in the years to come. 
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